Layering transitions for adsorbing polymers in poor solvents 
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An infinite hierarchy of layering transitions exists for model polymers in solution under poor 
solvent or low temperatures and near an attractive surface. A flat histogram stochastic growth 
algorithm known as FlatPERM has been used on a self- and surface interacting self-avoiding walk 
model for lengths up to 256. The associated phases exist as stable equilibria for large though not 
infinite length polymers and break the conjectured Surface Attached Globule phase into a series of 
phases where a polymer exists in specified layer close to a surface. We provide a scaling theory for 
these phases and the first-order transitions between them. 
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With the advent of sophisticated experimental tech- 
niques m, such as optical tweezers, to probe the be- 
haviour of single polymer molecules and the explosion 
of interest in the physics of biomolecules such as DNA 
there is a new focus on the study of dilute solutions of 
long chain molecules. It is therefore appropriate to ask 
whether the thermodynamic behaviour of such long chain 
molecules is well understood over a wide range of solvent 
types, temperatures and surface conditions. Even if one 
considers a fairly simple lattice model consisting of a self- 
avoiding walk on a cubic lattice with nearest-neighbour 
self-interactions in a half-space with the addition of sur- 
face attraction, the phase diagram has not been fully 
explored. In this letter we examine the whole phase dia- 
gram highlighting a surprising new phenomenon. In par- 
ticular, we demonstrate new features for large but finite 
polymer lengths involving the existence of a series of lay- 
ering transitions at low temperatures. 

Separately, the self-attraction of different parts of the 
same polymer and the attraction to a surface mediate 
the two most fundamental phase transitions in the study 
of isolated polymers in solution: collapse and adsorp- 
tion. Without a surface, an isolated polymer under- 
goes a collapse or coil-globule transition Q from a high 
temperature (good solvent) state, where in the infinite 
length limit the polymer behaves as a fractal with di- 
mension df = 1/u where u ^ 0.5874(2) (known as the 
extended phase) to a low temperature state {collapsed) 
where the polymer behaves as a dense liquid drop (hence 
a three-dimensional globule). In between these states is 
the well-known ^-point. Alternately a polymer in the 
presence of a sticky wall will bind (adsorb) onto the sur- 
face as the temperature is lowered j3|. At high tem- 
peratures only a finite number of monomers lie in the 
surface (desorbed) regardless of length even if the poly- 
mer is tethered onto the surface, while at low tempera- 
tures a finite fraction of the monomers will be adsorbed 



onto the surface in the large length limit and the poly- 
mer behaves in a two-dimensional fashion (with a smaller 
fractal dimension of 4/3). Much theoretical and experi- 
mental work has gone into elucidating these transitions. 
The situation when both of these effects are at work 
simultaneously and hence compete has received atten- 
tion in the past decade 0? SIE 0- three dimen- 
sions, four phases were initially proposed: Desorbed- 
Extended (DE), Desorbed-Collapsed (DC), Adsorbed- 
Extended (AE) and Adsorbed-Collapsed (AC) phases. In 
the AC phase the polymer is absorbed onto the surface 
and behaves as a two-dimensional liquid drop. Recently 
a new low-temperature (surface) phase named Surface- 
Attached Globule (SAG) " has been conjectured from 
short exact enumeration studies and the analysis of di- 
rected walk models ^3 • ^^i^ phase the polymer would 
behave as a three-dimensional globule but stay relatively 
close to the surface. In fact the claim is that there is not 
a bulk phase transition between DC and SAG (if SAG 
exists) in that the free energy of SAG and DC are the 
same. However, the number of surface monomers would 
scale as n^/^, where n is the number of monomers, rather 
than the as normally occurs in the desorbed state |J] . 

To explore the phase diagram it is natural to conduct 
Monte Carlo simulations. However the scale of the en- 
deavour becomes clear for even small system sizes (poly- 
mer lengths) because one is required to scan the entire 
two-energy parameter space of self-attraction and sur- 
face attraction. Even if once fixes one of the parameters, 
the study of the properties of the model as the other 
parameter is varied usually requires many simulations. 
Fortunately a recently developed algorithm ^Jl? Flat- 
PERM, is able to collect the necessary data in a single 
simulation. The power of this approach cannot be un- 
derestimated. We have utilised FlatPERM to simulate 
a self-avoiding walk model of a polymer with both self- 
attraction and surface attraction that allows us to calcu- 
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late quantities of interest at essentially any values of the 
energies. This was done with one very long simulation 
run for polymer lengths up to length Umax = 128 and 
also in multiple shorter CPU time runs, up to polymer 
length Umax = 256. We used a Beowulf-cluster to run 
these multiple simulations simultaneously with different 
random number 'seeds'. This allowed us to gain some 
estimate of statistical errors. 

The model 5j considered is a self-avoiding walk in 
a three-dimensional cubic lattice in a half-space inter- 
acting via a nearest-nearest energy of attraction Sb per 
monomer-monomer contact. The self- avoiding walk is at- 
tached at one end to the boundary of the half-space with 
surface energy per monomer of Eg for visits to the inter- 
face. The total energy of a configuration (pn of length n 
is given by 

Eni^n) -mb{(Pn)Sb - ms{(Pn)Ss (1) 

and depends on the number of non-consecutive nearest- 
neighbour pairs (contacts) along the walk 7715 and 
the number of visits to the planar surface m^. For 
convenience, we define = Sb/kBT and Ps = 
Es/kBT for temperature T and Boltzmann constant 
kB- The partition function is given by Zn{/3i)^/3s) = 

density of states. It is this density of states that is es- 
timated directly by the FlatPERM simulation. Our al- 
gorithm grows a walk monomer-by-monomer starting on 
the surface. We obtain data for each value of n up to 
^maxi and all permissible values of and m^. The 
growth is chosen to produce approximately equal num- 
bers of samples for each tuple of (n, m5,ms). The equal 
number of samples is maintained by pruning and enrich- 
ment . For each configuration we have also calculated 
the average height above the surface. Instead of relying 
on the traditional specific heat we have instead calcu- 
lated, for a range of values of and /3s, the matrix of 
second derivatives of \og{Zn{Phi Ps)) with respect to /^^ 
and (3s and from that calculated the two eigenvalues of 
this matrix. This gives a clear picture of the phase di- 
agram and allows for the accurate determination of the 
multicritical point 0| that exists in the phase diagram 
(see Figure 1(a)). 

We begin our discussion by showing a plot obtained 
from one run of the FlatPERM algorithm for Umax = 128 
in the region of parameter space that has been inves- 
tigated in previous works, namely < /^^ < 1.4 and 
< /^s ^ 1-6 (see Figure 1(a)). The phase boundaries 
seen by Vrbova and Whittington 6] are clearly visible 
with four phases in existence. For small /^^ and ps the 
polymer is in a desorbed and expanded phase (DE). For 
larger jSg adsorption occurs into the AE phase while for 
larger /^^ a collapse transition occurs into a phase de- 
scribed as the DC by Vrbova and Whittington Jp and 
either SAG or DC by Singh, Giri and Kumar 8]. We 
see little evidence for two phases in this region but given 



that the SAG/DC phase boundary is not a bulk phase 
transition this is not totally surprising! 

(a) 




FIG. 1: Plots of the logarithm of the largest eigenvalue of the 
matrix of second derivatives of the free energy with respect 
to pb and Ps for a range of bulk and surface energies. The 
lighter the shade the larger the value. The first figure shows 
the range of energies previously considered, and a schematic 
phase diagram consistent with Vrbova and Whittington. The 
white circle denotes the multicritical point. The second shows 
an extended range which clearly shows a new phenomenon (see 
main text). These plots were produced from our n < 128 long 
simulation run, using the n = 128 length data. 

The power of the FlatPERM method is that it al- 
lows us to explore regions of parameter space usually 
unavailable to canonical approaches and so we can con- 
sider a much wider range of Pb and pg. Of course the 
price paid is that the polymer lengths attainable are re- 
stricted due to both computer memory required and time 
needed to produce the samples. In Figure 1(b) we con- 
sider 0< pb< 5.0 and < < 5.0. 

To understand what is going on let us consider the 
mean density of surface contacts {ms)/n {coverage) as a 
function of Ps at fixed P^ = 4.0 for various values of n up 
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to Umax = 256 (Figure EJ. For small Ps the coverage is a 
slowly varying function of Ps and stays that way as n in- 
creases. For /3s larger than approximately /3b the coverage 
converges to a plateau of 1. So for /3s > /3}y essentially all 
the monomers are in the surface and the polymer should 
behave in a two-dimensional fashion. The transition to 
the maximum coverage regime (fully adsorbed) is quite 
sharp and reflects a first order phase transition in the 
thermodynamic limit: we shall confirm this inference be- 
low. The new phenomenon concerns intermediate values 
of /3s where other plateaus form at around a coverage of 
1/2 and, for larger n, also at 1/3. The transition from one 
plateau to another moves towards as n increases and 
also becomes sharp. We can interpret these intermediate 
"phases" as situations where the polymer is distributed 
roughly equally amongst a number of layers. For exam- 
ple when the coverage is 1/2 the polymer exists equally 
on the surface and in the layer one unit above the sur- 
face. As n increases more and more layer phases appear 
where the polymer exists in the first £ layers above the 
substrate. 
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FIG. 2: A plot of the mean density of visits {ms)/n versus [3s 
at f3b = 4.0 for lengths 64, 91, 128, 181, 256 (left to right) with 
statistical error estimated as shown. 

To confirm this picture let us consider the mean height 
of monomers above the surface (h) in Figure [31 for the 
same value of /^^ = 4.0. Assuming a uniform density 
across layers, the mean number of layers that the poly- 
mer subtends, {£) = 2(/i) + 1, can be deduced. We have 
also estimated the end-point position and the maximum 
height of the polymer, and the data agrees very well with 
this assumption. The average height can be seen to de- 
crease as Ps is increased in a series of plateaux corre- 
sponding to the plateaux of coverage. For the range of 
Ps where the coverage is approximately 1/2 the average 
height is almost exactly 0.5 and the maximum height of 
a monomer is 1 (two layers). Hence the average number 
of layers is 2, just as our hypothesis predicts. 

To explain the phenomenon of the layering seen above 
let us examine the zero-temperature situation. For pos- 
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FIG. 3: A plot of the average height of the polymer per 
monomer {h)/n (right axis) versus Ps at pb — 4.0 for lengths 
64,91,128,181,256 (left to right) with statistical error esti- 
mated as shown. On the left vertical axis the corresponding 
average layer number {£) (assuming uniform density). 



itive self-attraction and surface attraction the polymer 
will take on some compact configuration touching the 
boundary. Consider a Hamiltonian (fully compact) con- 
figuration of fixed height £ tethered to the surface. In 
particular, consider a rectangular parallelepiped with 
square cross-section parallel to the surface of side-length 
w. Hence we have n = £w'^ and the total energy E£ (ig- 
noring contributions from edges and corners) for a £-layer 
configuration is 

Eeisb^Ss) - -2ebn ^ {eb - Ss)^ ^ 2ebVhl (2) 

The energy can be minimised for fixed n when 

£3/2 = (^1 - (3) 

Since the system can only have integer values of ^, a par- 
ticular integer value of ^ > 2 will be stable for a range of 
Ss of size 0(1/ y/n). The AC phase, using this argument, 
which is given by ^ = 1, is stable for Ss > Sb and for 
some values of £s < Sb given by the relation Q. As £s is 
increased at fixed Sb the system's energy is minimised by 
smaller values of ^. At a fixed value of £s the difference 
between the energies of (£ + l)-layers and ^-layers scale as 
{£s—Sb)n£~^{£^l)~^. It can be argued that non-uniform 
layers are not stable (consider the total surface area of 
a block of smaller width on top of an ^-layer system) so 
that the system jumps from {£ + l)-layers to ^-layers at 
some value of e^. Hence we deduce that when the sys- 
tem swaps from {£ + l)-layers to ^-layers there will be a 
jump in the internal energy. We expect that this will be 
rounded by entropic effects at finite temperatures. 

The relation based on the zero-temperature en- 
ergy argument predicts that the transitions coalesce at 
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Ps = /35 as n tends to oo. For finite temperatures tiie po- 
sition of tiie transition need not be exactly P^. Let us de- 
note tiie infinite n limit transition as occurring at /3s = /3g 
The thermodynamic limit will realise a sharp bulk first- 
order phase transition at /3s = /3^. At finite polymer 
lengths each of the layering transitions are rounded ver- 
sions of the zero-temperature jumps in the internal en- 
ergy. Hence the layering transitions should be rounded 
first-order type transition with specific heat that grows 
linearly with system size and a transition width that is 
0{l/n). Note that since the layer phases are stable for 
segments of the /3s line of the order of 1/ and the tran- 
sitions take place in a region of /3s of the order of 1/n, 
the transitions become sharper as n increases. Finally, 
if we consider sitting at fixed value of /3s close to /^^ and 
increase the polymer length, we should see a set of layer- 
ing transitions between phases of layer £ and £ -\- 1 with 

We have tested the conclusions of the above argument. 
At Pb = 4.0 we have estimated the position of the ther- 
modynamic limit transition to the AC phase (ie 1-layer 
phase) to be (3^ ~ 4.07 by extrapolating the peaks of 
the fluctuations in the number of surface contacts for the 
strongest transition (ie from 2-layers to 1-layer) against 
Xj y/n. Figure^ shows a scaling plot of the logarithm of 
the fluctuations per monomer divided by n in the number 
of surface contacts against {P^ — Ps)V^' (For convenience 
we use the logarithm to display the needed scale.) The 
range includes the peaks from the 1-layer to 2-layers tran- 
sition, the 2-layers to 3-layers transition and the 3-layers 
to 4-layers transition. This demonstrates the scaling col- 
lapse of the height and shift of the layering transitions 
from ^-layers to ^+ 1-layers as ps is decreased. The shifts 
of the peaks of all three transitions scale towards the same 
estimate of P^ when using this same scale, l/y/n. The 
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FIG. 4: A plot of the logarithm of the fluctuations per 
monomer divided by n in the number of surface contacts at 
Pb — 4.0 with the horizontal axis scaled as (ps ~ Ps)y/n- We 
have used Ps = 4.07. Shown are lengths 128, 181, 256. 

width of the transitions can also be shown to scale with 
1/n, re-enforcing the hypothesis of first order transitions. 



In this letter we demonstrate that the fundamen- 
tal model of collapsing and adsorbing polymers in 
three dimensions contains a new phenomenon at low- 
temperatures; at finite polymer lengths a series of 
(rounded) layering transitions exist. These transitions 
increase in number and become sharper as the polymer 
length increases. We note that while this model is a 
lattice model, low temperature layering transitions have 
been seen in off-lattice models [l3| and arise due to the 
types of compact configurations that can occur in the ide- 
alised or physical polymer. The ability to coat a surface 
with a fixed thickness of polymer may have experimental 
and technological applications. We provide a theoretical 
framework based on zero-temperature energy arguments 
which explain these transitions. The arguments predict 
that the transitions coalesce in the infinite length limit 
to leave a transition between a collapsed, but not macro- 
scopically adsorbed, polymer and a collapsed polymer 
which is fully adsorbed. 
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